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ABSTRACT 

Using theoretical formalism of Chung-In Um et.al. the radial 

distribution function g(r) for two-dimensional liquid He 
interacting with Aziz potential has been evaluated for various 
densities. Our theoretical result indicates that g(r) increases with r 
and attains some maximum value and then becomes almost 
constant. 



INTRODUCTION 

1,2 

In an earlier paper we have 
evaluated the radial distribution function g(r) 
and ground state energy of two-dimensional 

liquid ^He interacting with Aziz potential 
and Leannard-Jones potentials using theor- 

3 

etical formalism of Chung-In Um et. al. In 

3 

this paper, using the same formalism , We 
have evaluated the radial distribution 

function of g(r) of two dimensional J He 
interacting through the Aziz potential at 
various densities. 

For many years, physicist have used 
hypernetted chain (HNC) and Fermi 
hypernetted chain (FHNC) theories to 
investigate the properties of highly 
correlated interacting boson and fermion 

systems such as liquid ^He and normal 



liquid 3 He. 4 ' 5 The HNC/FHNC equations 
with Jastrow ground-state wave functions 
treat consistently both long and short 
correlations of the system and give 
agreements to some degree with 
experimental results. In spite of qualitative 
successes of the Jastrow ground-state 
function, there exist quantitative differences 
between HNC/FHNC results and 
experimental ones, especially the 
equilibrium ground-state energy, and this 
fact leads one to consider the contributions 
arising from the three-body correlation 
functions. It is well known that the three- 
body correlations do not affect seriously the 
structural properties of ground states, such 
as the radial distribution function and liquid 
structure function, but contribute 
significantly to the ground-state energy. 

Since liquid J He system is less dense than 
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liquid He, that is, the equilibrium density of 
the former is 0.0166 A while that of the 

o -3 

latter is 0.02185 A in a three-dimensional 
system, effects of three-body correlations in 
3 He are rather small compared to those in 

4 

He. At equilibrium density, they make up 
about half the difference between the 
Jastrow results and experimental energy in a 
liquid 3 He system and lower the Jastrow 
ground-state energy by about 10 % in liquid 

3 He.' 

A trial ground-state wave function 
including three-body correlation functions 
in addition to two-body correlation effects 
can be constructed in the form 

V (r 1 ,r 2 ,...f N ) = n/4r lj )n/ 3 (r ij ,r jk ,r kl )|*) = F|<I.) 

i<i i<j<k 

where f 2 (ry) is the two-body correlation 
function which describes spatial correlations 
and depends only upon the spatial distance 
ry between two particle f 3 (ry, r jk , r Ki ) is a 

three-body correlation function, and I <I>) is 
the anti-symmetric product of non- 
interacting single-particle plane waves 
normalized in the system's area Q. with N 
fermions, which becomes simply unity for 
boson systems. 

In addition to the three-body 
correlation functions in the ground-state 
wave function, one should take into account 
the contributions coming from the 
elementary diagrams as well as effects of the 
HNC/FHNC equations to obtain better 
results. Unfortunately, there does not exist a 
general closed expression for the elementary 
diagrams and therefore one must calculate 
each diagram individually, which is very 



time consuming. During the last two 
decades, two approximations for the 
estimation of the elementary diagrams have 
been developed. One is known as the scaling 

6 

approximation, and the other is the 

7 

interpolating equation approximation. Both 
techniques have shown their efficiency for 
central two-body correlations, providing 
similar results for the total ground-state 
energy to those of the variational Monte 
Carlo simulation. The scaling approximation 
uses the fact that successive contributions of 
the five, six and higher order elementary 
diagrams are approximately proportional to 
and have very similar spatial behaviour to 

6 

the four-body elementary diagrams, so that 
one can represent the total elementary 
diagram contribution by means of the readily 
calculable four-body elementary diagrams. 
The interpolating equation approximation 
determines a parameter such that it gives the 
same value for the classical isothermal 
compressibility obtained from both the 

HNC/FHNC equations and Percus-Yevick 8 
equations using the classical pressure 
derivative or the compressibility integral. In 
this paper, one uses the scaling 
approximation to calculate the contribution 
of the elementary diagrams. 

Many works have been reported 
regarding the properties of three- 
dimensional ^He and ^He systems using the 
variational HNC/FHNC methods, the 

variational Monte Carlo method and the 
Green's function Monte Carlo method 

9 

(GFMC), there are few papers on two- 
dimensional systems. Using the Jastrow 
wavefunctions and solving an Euler- 
Lagrane equation within the HNC 
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10 

approximation. Hatzikonstinou studied the 

ground-state of two-dimensional liquid ^He 

interacting through a Lennard-Jones 

potential and explained the long-range and 

short-range behaviours of the correlation 

function to a certain extent. However, since 

he ignored the effects of the three-body 

correlation and elementary diagrams. His 

results give ground state energies too high 

n 

compared to the GFMC and diffusion 

12 

Monte Carlo (DMC) results. Chang 
calculated the ground-state energy and 

structure functions of two-dimensional ^He 

13 

in liquid density ranges by taking the same 
potential as an interatomic potential and 
using a self-consistent paired-phonon 
analysis. He included the effects of the 
three-body correlation in the ground-state 
wavefunction, but did not consider the 
elementary diagrams at all. Therefore, 
though he obtained very optimized liquid 
structure functions and radial distribution 
functions, his results for the ground state 
energy are still higher and the equilibrium 
density is lower than those of the GFMC and 
DMC calculations. 

For the ^He system, Novaco and 

14 

Campbell obtained theoretically the 
ground-state energy of quasi-two- 
dimensional helium absorbed on graphite by 
introducing a trial wavefunction which has 
only two-body correlations, finite extent 
orthogonal to the substrate, and the same 
translational symmetry as the substrate. 
Furthermore, they considered the effects of 
Fermi Dirac statistics to only the lowest- 
order correction, which is the energy of the 
Fermi sea. Using a variational calculation 
and the quantum theorem of corresponding 



states, Miller and Nosanow analysed 
properties of the two-dimensional ^He 

15 16 

system. Recently Brami et. al. investi- 
gated two-dimensional ^He system 
absorbed on graphite through variational 
Monte Carlo calculations considering a one- 
body wavefunction to describe a z- 
delocalization of helium atoms in addition to 
planar correlations in the total wavefunction 
which has only two-body terms. 

MATHEMATICAL FORMULAE USED 
IN THE EVALUATION 

The Hamiltonian for the considered 
system is given by 

H = T + V = XT i+ £v; j =-AXV? + 2v(|r 1 -f 1 |) 

i i<j 111 i i<j 

(1) 

where V(r) is the interatomic interaction. A 
popular potential used in theoretical 
investigations is that of Lennard-Jones, 
which has two parameters, i.e., the hard- 
core radius and the well depth, and can be 
written as 





12 











6 



£ = 10.22 K, = 2.556 A (2) 

However, V^j(r) includes only the 

dipole-dipole interaction and does not take 
into account multiple interactions. A more 
physically realistic potential which accounts 
for the self-consistent field Hartree-Fock 
repulsion and multiple interactions is the 
HFDHE2 potential of Aziz et. al. ^ Presented 
as 
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V u (r) = e ' 



A exp 



-a 



+ c, 



F(r) 



(3) 



where 



F(r) = 



_ I exp 



r/r m <D 
r/r. <D 



e*/k B = 10.8 
C 6 = 1.3732412 
C 10 = 0.178100 
r m = 2.9673 A 



(4) 



-(Dr m /r-l) 
1 

The values of the constants are 

A = 0.54485 x 10 6 , 

a = 13.353384, 

C 8 = 0.4253785, 

D = 1.241314, 

It is well known from the Green's 
function Monte Carlo simulations and other 
variational calculations that the Aziz 
potential gives closer results to experiments 
than Vj^j in three dimensions. One uses both 

potentials and compare the results with those 
in two dimensions. 

The n-particle distribution function is 
represented as 

p n =(?„...?„) = v n N(N-l)...(N-l) 
J \|/ + (f lt ... f n )\|/(f 1 ,...f n )dr n+1 ...df N 



x 



(¥'¥> 



(4) 

where ?! denotes the spatial coordinates of 
the i-th particle and v is the spin degeneracy 
of the system (2 for this system). In the 
homogeneous system, as N — > °° and £1— >°°, 
the single-particle distribution function 
reduces to the density of the system so that 
the radial distribution function can be 
expressed interms of the density and two- 
particle distribution function P2O") as 



V 



g(r) = - 5 -p a (r) = N(N-l) 

r 



J\|/ + (r 1 ,...fJ\|/(f 1 ,...fJdf 3 ...df N 

(v'v) 



(5) 



One uses the FHNC approximation 
to sum the diagrams arising from the cluster 
properties of g(r). In this scheme, g(r) can be 
decomposed as 



g(r y ) = g dd (r y ) + 2g de (r y ) + g ee (r y ) 



(6) 



where dd, de and ee represent terms in 
which both i and j are not exchanged, only j 
is exchanged, and both i and j are 
exchanged, respectively. The components of 
g(r-) are given by 

gddfa) = exp[u 2 (r y ) + N dd (r y ) + E dd (r y )] 
g de (r y ) = g dd (r y )[N de (r y ) + E de (r y )] ^ 

eJtil =^j^V /v +N Jti) + K( r 9 ) + M4 2 

gee(r y ) = g dd (r y )L(r y )/v 
where 

L(r y ) = -l(k F r y ) + v(N ee (r y ) + E ee (r y )) 

(8) 

l(x) = 2j 1 (x)/x (9) 

k F is the Fermi momentum of the system, 
and j^x) is the Bessel's function of the first 
kind of order 1. N mm represent sums of the 
nodal diagrams and E nm sums of the 
elementary diagrams. The equation g ee (rij) 
denotes the terms in which both i and j are 
exchanged in an incomplete exchange loop, 
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and N ee (ry) and E ee (ry) are sums of the nodal dressed with chains due to the three body 
diagrams and the elementary diagrams in correlation functions, represented as 



which i and j belong to the same permutation 
loop, respectively. 

The sums of nodal diagrams are 
given by the following integral equation: 

N dd( r ij) = r 2 [(gdd +gde "N dd " N de - 1)^,(6*, - 1)J 

+r 2 [(g dd -N dd -i) ik ,(g de ) 



c = r 

Wd L 3 



(16) 



c de = r 3 



(Sdd ^§de)i k '(Sdd) k j 
(gee _ gde) ik '(§dd) k j + ^(Sde)^ >(§de) k j 

(17) 



(10) 



N 



dc (r y ) = T 2 [(g dd + gclc - N M - N dc - l) ]k , (g dc - l) k . ] 



c„ = r, 



(2g de gde) ik >(gde) kj 
(§ee)i k '(gee)kj 



+r 2 [(g 



N 

dd x> dd 



l)ik'(8ee) k j 



(11) 



N ee(r„) = r 2 [(g de + g ce -N de -N ee ) ik ,(g de ) k . 



[(gde N de ). k ,(g ee ) ]i 



(12) 



N ee (r,) = r 2 [(g ee -N ee+ 1/V) ik ,(g ee ) kj ' 

(13) 

where r [ x ik>yki] is the convolution 
integral of x lk and y kj defined as 



(18) 

(19) 
where 

r 3 (x lk , y M ) = p \ d 2 r k [f 3 (r y , r jk , r B ) - l]x(r ik )y(r kj ) 

(20) 

For numerical calculations, one 
adopt the McMillan type function as a trial 
two body correlation function, which is 
widely used in variational and Monte Carlo 
calculations for liquid and solid helium 
systems : 



u 2 (r) = 



J (21) 

where b is a variational parameter chosen to 

r 2 ( x i k 'y k j) = P J d r k x(r ik )y(r ki ) ^4) give an optimised radial distribution 

function at a given density. To obtain more 

u 4.1. u ■ ii 11 i precise results, one include three body 
where the sum over each is over all allowed F , . „„ ' „ J 

combination of x, y, z, z', y' and z' which 

form proper exchange, and the sums of the 

Abe contributions in the FHNC 

approximation are neglected. The Abe where 

contributions will be taken into 

consideration when we include the 

elementary diagrams through the scaling 

approximation. To include three body 

correlation effects, one replace N nm (ry) + 



correlation effects in the following forms: 
u 3 (r y , r & , r u ) = £ T)(r y Jr^ )(?„ . f ft ) 



eye 



r((r) = V^i r ex P 



r ~ r i 



(22) 



(23) 



E nm (fij) in Eqs. (7) and (8) by N nm (r,j) + E nm 



where eye denotes the cyclic permutation 
among the three particle coordinates, i, j and 



(ry) + C mn (ry). Here C mn (ry) are diagrams k and 'J is a unit vector along the line 
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Table Tj 

An evaluated result of g(r) as a function of r(A) for two dimensional liquid 3 He interacting through Aziz 

Potential at various densities 



r(A) 


g(r) 


IVktPTltifll 

density 
0.01 A-2 


l IJldllltll 

density 
0.015 A-2 


l (JIVlllltll 

density 
0.02 A-2 


l (JIVlllltll 

density 
0.025 A-2 


2.0 


0.252 


0.307 


0.326 


0.355 


2.5 


0.296 


0.358 


0.389 


0.418 


3.0 


0.348 


0.406 


0.433 


0.453 


3.5 


0.409 


0.477 


0.507 


0.533 


4.0 


0.486 


0.532 


0.556 


0.588 


4.5 


0.585 


0.613 


0.639 


0.647 


5.0 


0.697 


0.722 


0.758 


0.768 


5.5 


0.786 


0.806 


0.844 


0.855 


6.0 


0.865 


0.952 


0.978 


0.998 


6.5 


0.954 


0.997 


01.067 


1.069 


7.0 


0.908 


1.052 


1.087 


1.092 


7.5 


0.885 


1.023 


1.042 


1.056 


8.0 


0.842 


1.007 


1.032 


1.035 


9.0 


0.807 


0.987 


1.007 


1.015 


10.0 


0.787 


0.973 


0.988 


1.002 



connecting particles i and j. Since it was first 
used in a variational Monte Carlo 
calculations, 12 this form for the three body 
correlations is generally used in variational 
HNC/FHNC methods. 13 The above 
parameters X u r { and C0i can be determined 
through a variational procedure also. In the 
three dimensional calculations, the value of 
the parameters from the Monte Carlo 
simulations are used, but one should 
determine these parameters for a two 
dimensional system through a variational 
procedure and HNC scheme because they 
are different in the two systems and the 



Monte Carlo results do not exist in two 
dimensions. They depend very weakly on 
the density, so that one use values at the 
equilibrium density for all density ranges. 

DISCUSSION OF RESULTS 

In this paper, we have evaluated the 
radial distribution function g(r) as a function 
of r for two dimensional liquid 3He 
interacting through Aziz potential at various 
densities. The evaluation has been 
performed with the help of Chung-In-Um 
et.al. 3 formalism. Our theoretical calculation 
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indicates that g(r) increases with r and 
attainis some maximum value at some value 
of r and thereafter it becomes almost a 
constant value. This behaviour is something 
different from g(r) of liquid 4 He in which 
g(r) increases and attains some maximum 
value and then decreases. 
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